Decay of Rabi oscillations by dipolar-coupled dynamical spin environments 
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We investigate the Rabi oscillations decay of a spin decohered by a spin bath whose internal 
dynamics is caused by dipolar coupling between the bath spins. The decay form and rate as a 
function of the intra-bath coupling is studied analytically, and confirmed numerically. The decay 
in general has neither exponential/Gaussian or power- law form, and changes non-monotonically 
with the intra-bath coupling, decelerating for both slow and fast baths. The form and rate of Rabi 
oscillations decay can be used to experimentally determine the intra-bath coupling strength for a 
broad class of solid-state systems. 
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Measurement of the Rabi oscillations decay is an im- 
portant step in studying decoherence of quantum sys- 
tems. For instance, extensive studies of Rabi oscillations 
in superconducting qubits P, 0, Q greatly enhanced un- 
derstanding of the decoherence caused by the bosonic 
baths and by the 1// noise. The most promising di- 
rections for improvement were identified, and the deco- 
herence time for superconducting qubits have been ex- 
tended to microseconds range placing them among 
most promising solid-state qubit candidates. Recently, 
much progress has been achieved in experimental im- 
plementation of long-living coherent Rabi oscillations in 
various spin systems: magnetic molecules [5|, NV impu- 
rity centers in diamond [y, 0i Q j rare-earth ions in non- 
magnetic host 9] , quantum dots 13, III, [l2| , to mention 
a few. Many of these systems are very attractive for 
basic studies of quantum spin coherence effects, and also 
constitute promising candidates for quantum information 
processing, coherent spintronics, or high-precision mag- 
netometry devices, provided that detailed understanding 
of the decoherence processes will be achieved. 

A major decoherence source in these spin systems is 
the coupling of the central spin (e.g. the electron spin of 
the NV center in diamond) to other spins present in the 
sample (environmental bath spins, e.g. the spins of ni- 
trogen atoms in diamond). Moreover, for many relevant 
spin systems, the direct coupling between the environ- 
mental spins is essential, producing internal dynamics of 
the bath. Here we explore theoretically the influence of 
such a dynamical spin bath on the Rabi oscillations of the 
central spin. Avoiding the commonl y u sed framework of 
generalized Bloch equations l^, 13,13, [Hi , we are able 
to investigate the form and rate of decay as a function 
of the intra-bath coupling strength. We find interesting 
behavior of the Rabi oscillations, which contradicts the 



expectations based on standard Redfield-type analysis: 
e.g., the slow bath leads to pronounced decay, while for 
the fast bath the decay rate vanishes but the Rabi fre- 
quency becomes renormalized. We demonstrate how the 
form and rate of the Rabi oscillations decay can be used 
in experiment to characterize the intra-bath dynamics, 
and provide a rather simple recipe for analysis of exper- 
imental data. 

Specifically, we consider the situation of dilute dipolar- 
coupled bath: in a non-magnetic host crystal, a single 
central spin of species S (e.g. belonging to a paramagnetic 
impurity) is coupled to a bath of dilute spins of species / 
(e.g. belonging to other kind of paramagnetic impurities 
or to nuclear spins) . The coupling between the bath spins 
is non-negligible, and is caused by dipolar interactions. 
This situation encompasses a wide range of interesting 
solid-state spin systems, from Er ions in CaWOiStudied 
in 1960s [i3] to the NV centers in diamond 0,0, Si which 
gained much attention very recently. Moreover, experi- 
mental progress now makes possible detailed studies of 
Rabi oscillations of an individual S spin [3, [Til . We 
assume that a large magnetic field is applied along 
the z-axis (as in standard NMR/ESR settings), leading 
to Zeeman splittings ujs = Js^z for the central spin S 
and LUj — "fiBz for the bath spins Ik (75 and 7/ are the 
gyromagnetic ratios of the S and / spins, respectively). 
Also, strong Rabi driving field Hji is applied at the fre- 
quency ids- The difference \lus — w/| is much larger than 
any other energy scale, and after transformation into ro- 
tating frame we obtain a standard secular Hamiltonian 
for two dipolar-coupled spin species [H, : 

H = h,S^ + AkSzPk + /f - lull) (1) 
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where 5"^'!'- 



and 1^'^'^ 



are the spin operators in the ro- 
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tating frame, and = Hji/2 is the rotating-frame Rabi 
driving field. The coupling constants Ak = 7s7/[l — 
3(n^)^]/r^ are determined by the positions r^. of the 
bath spins Ik {k = 1,...N), where = |rfc| and 
= J^k/ fk (the origin of the coordinate frame coincides 
with the central spin). Similarly, the intra-bath couplings 
Cki = 7/[l ^ 3(n^;)^]/r^.; are determined by the vectors 
Tfci — Tfc ^ I";. Note that the same Hamiltonian ^ can 
be obtained without external field Bz if the transition 
frequencies us and w/ are determined by the zero-field 
splitting (e.g. due to anisotropic interactions). Similarly, 
the assumption (used below) that S and Ik are spins 1/2 
is not essential: for larger spins, we can consider each 



pair of levels as a pseudo-spin 1/2 lie 

Initially the central spin is in the state " up" , and 
the bath is in a maximally mixed state (unpolarized 
bath at high temperatures), i.e. the initial density 
matrix of the whole system is p{0) = 2^^\ t)(t | (g) 



1b, where Is is the 2 x 2 identity matrix. This 
is appropriate for most experiments (for nuclear spins 
at temperatures above few nK, for electron spins — 
above tens of K). We calculate the time-dependent z- 
component of the central spin {Sz{t)) = TTp{t)Sz, where 
p{t) — exp {~iHt)p{0) exp (iHt). In order to see well- 
pronounced long-living Rabi oscillations, the driving field 
should be large, so we assume that hx is much larger 
than all other energy scales. We calculate the oscillations 
damping in the lowest order in 1/hx, treating the second 
(spin-bath coupling) and the third (bath internal Hamil- 
tonian) terms in Eq. [T] perturbatively, and excluding the 
bath internal Hamiltonian Hb ~ X]fe i ^ki^Iklf — Ifclz) 
by the interaction representation transformation. The 
resulting second-order Hamiltonian 



Bn = KSx + S^,B'(£)l(2K) 



(2) 



where B(t) = exp {iBBt)B exp {—iBBt), and the opera- 
tor B = j:kAkII. 

The evolution of B{t) is complex, involving intricate 
correlations between bath spins. However, exact dy- 
namics of every single bath spin is not important, since 
{Sz{t)) involves tracing over all bath spins. This is typi- 
cal for many spin-bath decoherence problems, from mag- 
netic resonance to quantum information processing, and 
many approaches have been dev elop ed from the early 
days of NMR/ESR theory [11, ^ till very recently 
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25|, |26|, |27|, |28l . All approaches involve 
a trade-off between quantitative rigour and qualitative 
understanding. Below, following the works [H, [l^l, we 
approximate the effect of the bath by a random field 
B{t), which is modeled as an Ornstein-Uhlenbek process 
with the correlation function {B{t)B{0)) = exp {—Rt), 
where the dispersion b = \/'^k ^fc' while the correlation 
decay rate R is determined by Bb- This model may be 
oversimplified for complicated situations, e.g. the descrip- 
tion of advanced contr ol p rotocols requires more sophis- 
ticated treatment 1241 12811. However, our direct numeri- 



cal simulations evidence that this model is quantitatively 
adequate for description of the Rabi oscillations decay, 
while providing a clear description of the physics under- 
lying the Rabi oscillations decay, and allowing access to 
the regimes outside of the Bloch equations framework. 

The Hamiltonian @ reflects a simple physical picture. 
The zero-order eigenstates of the Hamiltonian ^ corre- 
spond to the central spin states |-|-) = l/-\/2[| t) + I i)] 
and |— ) = 1/V2[| t) ^1 i)], separated by a large Rabi 
frequency hx- Since hx ^ b and hx ^ R, the field B{t) 
has no spectral components of noticeable magnitude at 
the Rabi frequency, and does not lead to transition be- 
tween the states |-|-) and |— ). The only relevant process is 
the pure dephasing, when the field B destroys the initial 
phase relation between the states |-|-) and |— ), leading to 
decay of {Sz{t)), so that 

(Szit)) ^ 1/2 (cos$) = 1/2 Re (exp (i$)) (3) 

where (•) denotes average over all possible realizations of 
B{t), and the phase 



1 /■* 

<^ = hxt+ — / B'^{s)ds = hxt + Q 
2hx Jo 



(4) 



is the total phase difference between the states |-|-) and 
|— ) accumulated during time t, cf. the Hamiltonian ([2]). 
The key quantity M{t) = (exp (i9)) is an analytically 
computable Gaussian path integral over the Ornstein- 
Uhlenbeck process, which gives the answer 



mt)]-' 



1/2 Re[M{t)exp {ihxt)] 
exp i-Rt) [cosh Pi + (R/P) sinhPi] 
. b^ 



hxP 



exp {—Rt) sinh Pt, 



(5) 



where P = ^ R^ - 2ib'^R/hx. 

Analysis of Rabi oscillations is often based on Bloch- 
type equations [3] or its generalizations EBl , derived for 
various systems from quantum optics [16[ to solid state 
fist - It is based on the Redfield-type approach 1^ 14 1, 
taking into account only the terms which are secular with 
respect to the Rabi driving hxSx- In addition to dephas- 
ing, these terms describe the actual transitions between 
the states |-|-) and |— ), which lead to a longitudinal re- 
laxation (along the x-axis) of the central spin with the 
rate F/ ~ b^R/hx- This rate is of second order in 1/hx, 
and is determined by the spectral density of B{t) at the 
Rabi frequency hx- Also, the generahzed Bloch equa- 
tions, having constant coefficients, always predict the de- 
cay to have a (multi)exponential form. In contrast, our 
results are not limited to the terms secular with respect 
to Rabi driving, and give the decay rate of the first order 
in l/hx- The solution ([5|) predicts the decay which has, 
in general, no simple form (multiexponential, Gaussian, 
power-law, etc.) The familiar exponential decay occurs 
only at special values of 6, hx, and R, see below. 
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The effect of the bath internal dynamics may be impor- 
tant even for slow baths, with R <C iP' /h^- In this limit, 
P = by/-2iR/h^, see Eq. [5l At short times t < l/\P\, 
the bath behaves as static, and the Rabi oscillations en- 
velope exhibits non-exponential slow decay of the form 

[l -I- {}Pt/hx)^\ "^'^^j in accordance with the exact results 
obtained earlier 0, [H, [H, [l^l . At long ti mes t > 1/|P|, 
the decay has exponential form exp {—htyjR/Ahx), with 
the rate which decreases very slowly with decreasing R. 
As a result, even for very slow bath the effect of the 
bath dynamics remains noticeable. In experiments, the 
above-described behavior of the Rabi oscillations decay 
can be detected by varying (since the ratio of R to 
b'^ /hx determines how fast/slow the bath is), and makes 
it possible to estimate b and R. 

Another interesting feature of our results is that the de- 
cay of Rabi oscillations changes non-monotonously with 
R: it is fastest for R ~ b'^/hx, slowing down for both 
slow and fast baths. This is confirmed by direct numeri- 
cal simulations, see Fig. [T] and discussion below. 

The regime of fast bath, with R ^ b^ /hx, is even more 
interesting and unexpected. In this case 



Mit) = exp 



bH 



bH 



2hx AhlR 



(6) 



and the Rabi oscillations exhibit exponential decay with 
the rate Ft = b"^/ (Ah^R), which vanishes quickly for large 
R. However, it does not mean that the effect of the fast 
bath disappears. The bath still noticeably affects the 
central spin, shifting its Rabi frequency by b'^ /2hx_, and 
only the decaying part of the bath effect vanishes [291 [sot. 
This slowing of the Rabi oscillations decay is confirmed 
by direct numerical simulations (Fig. [5]) . 

On the other hand, for fast baths, the transitions be- 
tween the states |+) and |— ) become important. While 
the dephasing between these states just shifts the Rabi 
frequency, the transitions lead to longitudinal relaxation 
along the x-axis with the rate F; ~ b^R/h^. This im- 
plies [l^l exponential decay of Rabi oscillations with the 
rate Ti/2, which is comparable with the decay rate Ft 
caused by the pure dephasing. Above, for simplicity, we 
omitted the longitudinal relaxation (since it is of order 
but we can include it using the Redfield-type ap- 
proach: the answer ([5]) for Sz{t) just has to be multipHed 
by exp(— F;t/2). This explains how the system enters 
the generalized Bloch equation regime at i? S> &: the 
dephasing effect becomes negligible and the longitudinal 
relaxation becomes dominant. In experiment, this regime 
can be identified by comparing F; and Ft : they would be 
of the same order, changing as for all sufficiently 

large hx- 

To ensure that the physical picture above is adequate 
for real dipolar-coupled bath, we performed direct numer- 
ical simulations of the Rabi oscillations decay. We place 
the central spin and N bath spins randomly, with uniform 
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FIG. 1: (Color online). Exact numerical simulations for 
= 15 bath spins, (a) — Correlation function {B{0)B{t)} 
(normalized to 1 at t = 0) obtained from direct simulations 
(black), see also [s^j- Its fitting (green) determines the bath 
parameters &i, 62, and R for each value of the intra-bath 
coupling scale Eb (here Eb = 1, fei = 0.62, &2 — 0.58, 
R = 0.095). (b) — Rabi oscillations decay for = 14.14, 
b = 0.85, and Eb ~ 1; individual oscillations are not resolved. 
Numerical results (black) agree well with analytics (magenta 
line, only the envelope shown), (c) — Envelopes of simulated 
Rabi decay for Eb = 1 (blue), Eb =0.1 (green), and Eb — 
(static bath, red); corresponding individual oscillations near 
t = 300 are shown on panel (d) by same colors. Analyti- 
cal results practically coincide with simulations, and are not 
shown. The decay rate changes nonmonotonically with Eb'- 
the decay is slower for Eb = 1 (blue), and Eb = (red) in 
comparison with Eb ~ 0.1 (green). At t = 300, on panel (d), 
the oscillation amplitude for Eb = 0.1 (green) is twice smaller 
than that for Eb ~ 0. 



density n = 1, inside a cube with the side (iV-f 1)^/'^, cen- 
tral spin being in the center of the cube. All interaction 
coefficients are calculated according to Eq. [U using the 
actual coordinates of the spins. The value Esb = 7s7/j 
which determines the strength of coupling between the 
central spin and the bath, is set to 1, thus defining the 
energy and time scales for all quantities below. The value 
Eb = 7I, which governs the energy scale of intra-bath 
couplings, is varied, making the bath slower or faster. We 
simulate the dynamics of the system using two numeri- 
cal approaches. For small N (e.g., iV = 15 in Fig. [T]) we 
exactly solve of the Schrodinger equation with the Hamil- 
tonian ^ via Chebyshev polynomial expansion [sij . The 
dipolar-coupled systems with iV > 15 are difficult to 
model this way (due to exponentially increasing resources 
requirements), so we use the P-representation sampling 
[sit for modeling larger baths {N = 159 in Fig. To 
compare numerical solutions with the analytics, we cal- 
culate the total couphng to a bath b = EfcLi^fc]^^^ 
directly from the positions of the spins (see Eq. [1] and 
below) . The correlation decay rate R requires a separate 
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ity of the approximation is confirmed by direct numer- 
ical simulations. The oscillations decay has interesting 
features, such as non-monotonic variation of the decay 
rate with increasing the intra-bath coupling, and slow- 
ing down of the decay for fast baths. Studying the be- 
havior of the Rabi oscillations may help in experimen- 
tal characterization of the dynamical spin bath, and is 
well within experimental reach e.g. for NV centers in di- 
amond. Work at Ames Laboratory was supported by 
the Department of Energy — Basic Sciences under Con- 
tract No. DE-AC02-07CH11358. We acknowledge sup- 
port from AFOSR (D.D.A.), FOM and NWO (R.H.). 
A.E.F. acknowledges support from the Microsoft Corpo- 
ration. 



FIG. 2: (Color online). Numerical simulations for — 159 
bath spins, (a) — Correlation function {B{0)B{t)) (normal- 
ized to 1 at t = 0) obtained from P-representation simula- 
tions (black); its fitting (green) determines the parameter 
R for each value of the intra-bath coupling scale Eb (here 
Eb ~ 1 and R = 0.097). In spite of rather large statis- 
tical fluctuations, the parameter R is determined precisely 
enough to be used in Eq. [S] (b) — Rabi oscillations decay 
for hx = 15.0, b = 1.39, and Eb ~ 0.1; individual oscilla- 
tions are not resolved. Numerical results (black) agree well 
with analytics (magenta line, only the envelope shown), (c) 
— Envelopes of simulated Rabi decay for Eb =0.1 (blue), 
Eb ~ 1 (green), and Eb = 10 (red); the longitudinal decay 
for Eb = 10 is taken into account. Analytical results prac- 
tically coincide with simulations, and are not shown. The 
decay rate decreases for faster baths: the decay is slowest 
for Eb ~ 10 (red). Individual oscillations for Eb ~ 10 near 
t — 20 are shown on panel (d) in red: their frequency is 
shifted by /2hx — 0.064 from the value hx = 15.0 (oscilla- 
tions with frequency hx — 15.0 are shown by dotted black line 
to demonstrate the phase difference accumulated since t — 0). 



simulation: we calculate {B{0)B{t)) and find R by fitting 
it to a decaying exponent, see also [32|. Varying the sys- 
tem parameters in a wide range, we simulated baths with 
N = 15, 59, and 159 spins, and found good agreement 
between numerics and analytics; typical results are given 
in Figs, [m 

The resulting experimental recipe is rather simple. If 
the decay has a power-law form at shorter times, chang- 
ing to exponent later, with the decay constants changing 
as l/hx and the duration of the two regimes varying with 
hx, then the bath is slow. If the Rabi oscillations decay 
is exponential, with the rate changing as and the 

frequency shift varying as 1/hx, then the bath is fast. 
The fast bath can be made slow by strong increase in h^- 

Summarizing, we studied the decay of the Rabi os- 
cillations of the central spin interacting with a dipolar- 
coupled dynamical spin bath. Approximating the effect 
of the bath as a random field (Ornstein-Uhlenbeck pro- 
cess), we find analytically the form of the decay. Valid- 
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